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Abstract 

Using the low-energy effective field theory for magnons, we systemati- 
cally evaluate the partition function of the 0(3) ferromagnet up to three 
loops. Dyson, in his pioneering microscopic analysis of the Heisenberg 
model, showed that the spin-wave interaction starts manifesting itself in 
the low-temperature expansion of the spontaneous magnetization of an 
ideal ferromagnet only at order T 4 . Although several authors tried to go 
beyond Dyson's result, to the best of our knowledge, a fully systematic and 
rigorous investigation of higher order terms induced by the spin- wave inter- 
action, has never been achieved. As we demonstrate in the present paper, 
it is straightforward to evaluate the partition function of an ideal ferro- 
magnet beyond Dyson's analysis, using effective Lagrangian techniques. In 
particular, we show that the next-to-leading contribution to the sponta- 
neous magnetization resulting from the spin-wave interaction already sets 
in at order T 9//2 - in contrast to all claims that have appeared before in the 
literature. Remarkably, the corresponding coefficient is completely deter- 
mined by the leading-order effective Lagrangian and is thus independent 
of the anisotropics of the cubic lattice. We also consider even higher-order 
corrections and thereby solve - once and for all - the question of how 
the spin-wave interaction in an ideal ferromagnet manifests itself in the 
spontaneous magnetization beyond the Dyson term. 
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1 Introduction 



In a landmark paper on the description of ferromagnets at low temperatures [H, 
Bloch introduced the concept of spin waves and identified them as the relevant low- 
energy degrees of freedom. As an immediate application, he evaluated the leading 
coefficient in the low-temperature expansion of the spontaneous magnetization: this 
term, corresponding to noninteracting magnons, is of order T 3//2 . As is well-known, 
various authors subsequently tried to find the leading term in this series originating 
from the spin-wave interaction, ending up with conflicting results: corrections to 
Bloch's law both of order T 7//4 and T 2 were found The situation remained 

rather unclear until Dyson, in his pioneering analysis of the thermodynamic behavior 
of an ideal ferromagnet [6j , showed that the previous results were wrong altogether and 
that the spin-wave interaction in the spontaneous magnetization starts manifesting 
itself only at order T 4 . 

Dyson's motivation was the apparent contradiction between the various results 
published in the literature. He successfully solved this paradox by setting up a fairly 
complicated mathematical machinery - in his own words [6j: "The method of the 
present paper settled the disagreement by showing that both calculations were wrong". 
In fact, as Dyson states in Ref. [6], a third calculation existed that was also in con- 
tradiction with the other two. 

Within the last few decades several articles have appeared dealing with the struc- 
ture of the series for the spontaneous magnetization beyond the Dyson term. Various 
authors, using different methods, have given their account on what the temperature 
power of the next-to-leading order term due to the spin-wave interaction should be 
and how the general structure of the series beyond Dyson should look like. Not all of 
these results that have appeared in the literature over time, however, as we discuss 
in more detail later on, are consistent with one another. Our main motivation is thus 
reminiscent of Dyson's, namely, to determine which one of these calculations yields 
the correct low-temperature expansion for the spontaneous magnetization of an ideal 
ferromagnet. 

Due to its mathematical rigor, Dyson's calculation is not easy to understand and 
the perturbative scheme developed for the evaluation of the partition function is fairly 
complicated. Indeed, after Dyson's analysis, many authors tried to reproduce and red- 



erive his result with alternative methods in a more accessible manner |7H17||. Among 
these references we would like to point out the paper by Zittartz [11], upon which 
Dyson comments "Zittartz replaced my cumbersome mathematics by a simple and el- 



egant construction" 18]. Surprisingly, not all authors were able to confirm Dyson's 
result: in particular, a new interaction term in the spontaneous magnetization of or- 
der T 3 started to haunt the literature 19H29| - later on, this term was recognized as 
spurious and Dyson's series was confirmed to be correct 
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While all these studies were performed within the framework of microscopic or 
phenomenological theories based on the Heisenberg model, in the present work, we 
will follow another approach which has the virtue of being completely systematic 
and model-independent: the method of effective Lagrangians. Within the effective 
Lagrangian framework, the structure of the low-temperature expansion of the spon- 



taneous magnetization was analyzed in Ref. [36| up to order T and Dyson's series 
was reproduced in a straightforward manner. In the effective language, as we will see, 
this corresponds to including Feynman diagrams for the partition function up to two 
loops. The effective analysis also readily demonstrated that there is no interaction 
term of order T 3 in the low-temperature series of the spontaneous magnetization. 

In the present work we go beyond Dyson's analysis and explicitly calculate the 
effect of the spin-wave interaction beyond T 4 in the spontaneous magnetization of an 
ideal ferromagnet. To the best of our knowledge, this is the first time that the structure 
of this power series is given in a fully systematic and rigorous way. Going beyond 
Dyson's analysis then means that, in the effective Lagrangian framework, we have 
to consider Feynman diagrams up to three-loop order in the perturbative expansion 
of the partition function. As it turns out, in the spontaneous magnetization of an 
ideal ferromagnet, the next-to- leading interaction term already sets in at order T 9 / 2 
- remarkably, the corresponding coefficient is completely determined by the two low- 
energy coupling constants of the leading-order effective Lagrangian C 2 e ^. It does not 
involve any higher-order effective constants from where the anisotropies of the 
cubic lattice start showing up. 

Although several authors have also discussed the structure of temperature pow- 
ers beyond the T 4 -term, their conclusions are in contradiction with the systematic 
effective field-theory approach and therefore erroneous. In particular, to the best of 
our knowledge, none of the existing calculations ended up with an interaction term 
of order T 9//2 , which in fact represents the leading correction to Dyson's result. 

Within the effective Lagrangian framework, we then analyze the general structure 
of even higher-order corrections in the spontaneous magnetization originating from 
the spin-wave interaction, and show that these are of order T 5 , T 11 / 2 , T 6 , ... - again 
contradicting earlier calculations that have appeared in the literature. 

The effective Lagrangian method is based on an analysis of the symmetry proper- 
ties of the underlying theory, i.e., the Heisenberg model in our case, and can univer- 
sally be applied to systems with a spontaneously broken symmetry. It is formulated in 
terms of Goldstone boson fields which represent the dominant low-energy degrees of 
freedom. The effective Lagrangian method is very well established in particle physics, 
where the low-energy effective theory for quantum chromodynamics - chiral perturba- 



tion theory - has been constructed a long time ago [37|, |38] . There, we are dealing with 
a spontaneously broken chiral symmetry and the corresponding Goldstone bosons are 
the pseudoscalar mesons. Spontaneous symmetry breaking is also a common phe- 
nomenon in condensed matter physics and the effective Lagrangian method has in 
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fact been transferred to this domain in Ref. 39]: Magnons and phonons, e.g., are 



the Goldstone bosons resulting from a spontaneously broken spin rotation symme- 
try 0(3) — > 0(2) and a spontaneously broken translation symmetry, respectively. In 
particular, the leading-order effective Lagrangian for the 0(3) ferromagnet was con- 
structed in Ref. j39[ and the extension to higher orders in the derivative expansion 



was performed in Refs. [36J, |40 



The paper is organized as follows. Since the systematic effective Lagrangian 
method is still not very well known within the condensed matter community, in Sec. [2] 
we give a brief outline of the method, having in mind the ferromagnet as specific 
system. In Sec. 13.11 we briefly review the evaluation of the partition function of an 
ideal ferromagnet up to order T 5 . We then go beyond Dyson's analysis and extend 
the evaluation to order T 11//2 in Sec. 13.21 While the renormalization up to order T 5 is 
straightforward, the handling of ultraviolet divergences at order T 11 / 2 is more involved 
and is considered in detail in Sec. HI The low-temperature expansion of the partition 
function and various thermodynamic quantities is given in Sec. |5j Our main result - 
the low-temperature series for the spontaneous magnetization of an ideal ferromagnet 
beyond Dyson's analysis - is presented in Sec. |6j Here we also compare our results 
with the condensed matter literature. While our conclusions are presented in Sec. 
details on the numerical evaluation of a specific three-loop graph are discussed in 
Appendix |A] 

We would like to provide the interested reader with a list of publications that deal 
with applications of the effective Lagrangian method to condensed matter systems. 
Applications to systems exhibiting collective magnetic behavior include spin-wave 
scattering processes 4l|, spin- wave mediated nonreciprocal effects in antiferromagnets 
|42| . antiferromagnets at finite volume 43-46| and finite temperature 47, 48], spin 



waves in canted phases [49j and antiferromagnets in two dimensions doped with charge 



carriers 



50H55|. Further applications include phonons 56], SO (5) invariance and high 



T c -superconductivity [57J as well as supersolids [58j . Pedagogic introductions to the 



effective Lagrangian method may be found in Refs. [59H65 



In particular, we would like to point out that in a recent article on an analytically 



66 , the 



solvable microscopic model for a hole-doped ferromagnet in 1+1 dimensions 
correctness of the effective field theory approach was demonstrated by comparing the 
effective theory predictions with the microscopic calculation. Likewise, in a series 
of high- accuracy investigations of the antiferromagnetic spi n- qu antum Heisenberg 



model on a square lattice using the loop-cluster algorithm [67H70j]. the Monte Carlo 



data were confronted with the analytic predictions of magnon chiral perturbation the- 
ory and the low-energy constants were extracted with permille accuracy. All these 
tests unambiguously demonstrate that the effective Lagrangian approach provides a 
rigorous and systematic derivative expansion for both ferromagnetic an antiferromag- 
netic systems. 
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2 Systematic Low-Energy Effective Field Theory 
for Ferromagnetic Magnons 



The effective Lagrangian method is based on a symmetry analysis of the underlying 
system. In the present case we study ferromagnets, which are described by the Heisen- 
berg model. Nevertheless, the effective field-theory predictions are model-independent 
and universal, as they are valid for any system displaying the same symmetries as the 
Heisenberg ferromagnet. Microscopic details of the system are taken into account 
through a few low-energy coupling constants in the effective Lagrangian. Symmetry 
does not fix the actual numerical values of these couplings - in general, these have to 
be determined experimentally or in a numerical simulation of the underlying model. 
Symmetry, however, does unambiguously determine the derivative structure of the 
terms in the effective Lagrangian. 

The most important symmetry in the present case is the spontaneously broken 
spin rotation symmetry: Whereas the Heisenberg model, 

Hp = — J S m ■ S n , J = const. , (2.1) 

n.n. 

is invariant under global 0(3) spin rotations, the ground state of the ferromagnet 
(J > 0) is invariant under the subgroup 0(2) only. According to the nonrelativistic 
Goldstone theorem Tl T5|| , we then have one type of spin- wave excitation - or one 



magnon particle - in the low-energy spectrum of the ferromagnet which obeys a 
quadratic dispersion relation. 

The interaction between an external constant magnetic field H = (0, 0, H), H > 
and the spin degrees of freedom is taken into account through the Zeeman term. In 
the corresponding extension of the Heisenberg model, 

H = Ho-ti^2^ n -H, (2.2) 

n 

the magnetic field couples to the vector of the total spin. The above Hamiltonian, 
defined on a cubic lattice with purely isotropic exchange coupling between nearest 
neighbors, represents what Dyson called ideal ferromagnet. 

Apart from internal symmetries we also have to consider the various space-time 
symmetries. Compared to particle physics where we have Lorentz invariance, the sit- 
uation is more complicated in condensed matter physics, because the center of mass 
system represents a preferred frame of reference. Moreover, the crystal lattice singles 
out preferred directions, such that the effective Lagrangian need not be rotation- 
ally invariant. In the case of cubic geometry, however, it has been shown that the 
anisotropies of the lattice start manifesting themselves at higher orders of the deriva- 



tive expansion [45] - the leading-order effective Lagrangian is thus invariant under 
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space rotations. Moreover, as the effective analysis refers to large wavelengths, it does 
not resolve the microscopic structure of the crystal: the system appears homogeneous 
and the effective Lagrangian is also invariant under translations. 

The idea underlying the construction of effective Lagrangians is straightforward 



76| : One writes down the most general expression consistent with the space-time 
symmetries and the internal, spontaneously broken symmetry G of the underlying 
system in terms of Goldstone boson fields U a (x), a = 1, . . . , dim(G)-dim(H), where the 
group H refers to the symmetry group of the ground state. The effective Lagrangian 
then consists of a string of terms involving an increasing number of derivatives or, 
equivalently, amounts to an expansion in powers of the momentum. Furthermore, the 
effective Lagrangian method allows to systematically take into account interactions 
which explicitly break the symmetry G of the underlying model, provided that they 
can be treated as perturbations. In the present case we will include a weak external 
magnetic field H. 

For the 0(3) ferromagnet, the leading-order effective Lagrangian is of order p 2 and 
takes the form 

C lff = ^fuT + ^ HU3 - \f%UW . (2.3) 

The two real components of the magnon field, U a (a = 1,2) are the first two compo- 
nents of the three-dimensional unit vector U l = (U a , U 3 ), which transforms with the 
vector representation of the rotation group. While the structure of the above terms is 
unambiguously determined by the symmetries of the underlying theory, at this order, 
we have two a priori unknown low-energy constants: the spontaneous magnetization 
£ and the constant F. The above Lagrangian leads to a quadratic dispersion relation 

uj(k) = ^ + 0(\kC) , 7=^, (2-4) 

obeyed by ferromagnetic magnons. It is important to note that one temporal deriva- 
tive (do) is on the same footing as two spatial derivatives (d r d r ) - in the derivative 
expansion, two powers of momentum thus count as only one power of energy or tem- 
perature: k 2 cx u,T. 

Dyson evaluated the low-temperature expansion of the spontaneous magnetization 
up to terms of order T 4 or, equivalently, the partition function up to order T 5 . This 
then means that, in the effective Lagrangian framework, we have to consider the 
expansion of the partition function up to order p 10 . This calculation was performed 
in Ref. [36[. In the present work, we go one step further and consider the expansion 
beyond Dyson's analysis, taking into account diagrams of order p 11 . As it turns out, 
the corresponding contributions lead to a spin-wave interaction term of order T 9//2 in 
the spontaneous magnetization. 

The effective Lagrangian method provides us with a simultaneous expansion of 
physical quantities in powers of the momenta and of the external fields. The essential 
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Figure 1: Feynman graphs related to the low-temperature expansion of the partition 
function for a ferromagnet up to order p 10 in dimension d—3+1. The numbers attached 
to the vertices refer to the piece of the effective Lagrangian they come from. Vertices 
associated with the leading term C? e ^ are denoted by a dot. Note that ferromagnetic 
loops are suppressed by three momentum powers in d—3+1. 



point is that, to a given order in the low-energy expansion, only a finite number 
of effective coupling constants and only a finite number of graphs contribute. The 
leading terms stem from tree graphs, whereas loop graphs only manifest themselves 
at higher orders in the derivative expansion 37]]. So the question arises as to what 
order in the effective expansion we have to go - i.e., how many derivatives in the 
effective Lagrangian we have to include and how many loops we have to consider - if 
we want to evaluate the partition function of a ferromagnet up to order p 11 . 

While loops are suppressed by two momentum powers in a Lorentz-invariant frame- 
work, it was shown in Ref. [36[ that loop corrections involving ferromagnetic magnons 
are suppressed by three momentum powers 0. Up to order p 10 , as depicted in Fig. (TJ 
we thus have to consider graphs which involve two loops at most and have to take into 
account pieces of the effective Lagrangian involving up to six derivatives. At order 
p 11 , as depicted in Fig. [21 three- loop graphs start to show up. At the same time we 
also have two one-loop graphs which involve vertices from higher-order pieces of the 
effective Lagrangian: Diagram lid contains a vertex from J~-lff, while diagram lie 
contains insertions from and £e//- These five graphs represent the additional 
diagrams we have to evaluate when we go one step beyond Dyson's analysis. 



We now address the question regarding the explicit structure of the pieces C 



e//' 



: We are considering the case of four space-time dimensions. If one lowers the spatial dimension, 
loops are less suppressed: Loops for ferromagnetic magnons in d=2+l, e.g., are suppressed by two 
momentum powers, while loops for antiferromagnetic magnons in d=2+l are suppressed by one 
momentum power only. 
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Figure 2: Feynman graphs related to the low-temperature expansion of the partition 
function for a ferromagnet at order p u in dimension d=3+l. The numbers attached 
to the vertices refer to the piece of the effective Lagrangian they come from. Vertices 
associated with the leading term C 2 eff are denoted by a dot. 



J~-eff an d ^eff- Phst °f au > note that there are no contributions to the effective 
Lagrangian leading to odd momentum powers: The pieces £-ljf, ^e//; ■ ■ ■ necessarily 
involve terms with an odd number of space derivatives like 

c abc e rst d r U a d s U b d t U c , (2.5) 

which are excluded by parity - parity is a discrete symmetry of the underlying Heisen- 
berg model that has to be respected by the effective Lagrangian. 

The next-to- leading order Lagrangian is thus of order p 4 . It contains terms with 
two time derivatives, terms with one time and two space derivatives, and terms with 
four space derivatives. The time derivatives along with the magnetic field, however, 



can be eliminated with the equation of motion, such that L\r, takes the form 36 



eff 

C\ n = kidrU^rWf + kidrU^sU 1 ) 2 + ACT* ALT , (2.6) 

where A denotes the Laplace operator in three dimensions. The next-to-leading order 
effective Lagrangian hence involves the three effective coupling constants h, I2 and ^3. 

An inspection of the diagrams in Figs. [T] and [2] reveals that insertions from and 
J~-lff only appear in one-loop graphs: the only terms we need are thus quadratic in the 
magnon field. Eliminating again time derivatives and terms involving the magnetic 
field, the pieces relevant for our calculation are 

C 6 eff = cilTAW , C 8 eff = dxirtflT . (2.7) 
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We conclude this section with a remark concerning effects induced by the anisotropy 
of the lattice. Regarding the cubic lattice, we have mentioned that the anisotropies 
start manifesting themselves at the four- derivative level: the pieces C*ff, C- & e ff and 
indeed contain additional terms - not displayed in Eqs.( l2.6p and (12. 7p - which 
are not invariant under space rotations, but still invariant under the discrete rotation 
and reflexion symmetries of the cubic lattice, such as 

d s d s u l d s d s u l . (2.8) 

s=l,2,3 

In the present analysis, however, we neglect these extra terms and assume space 
rotation invariance up to order p s . The conclusions of the present paper regarding the 
manifestation of the spin-wave interaction in the partition function are not affected 
by this idealization: According to Fig. El the interaction contribution beyond Dyson 
is determined by the three- loop graphs of order p 11 : These graphs only involve the 
leading-order Lagrangian C 2 eff which is perfectly invariant under space rotations. 



3 Evaluation of the Partition Function 



The low-temperature expansion of the partition function for the 0(3) ferromagnet 
was evaluated in Ref. [36( up to order p 10 . In Sec. 13. II we briefly review some essential 
features of that calculation. In Sec. I3.2l we then extend the evaluation of the partition 
function to order p 11 . For a review of the effective Lagrangian method at non-zero 
temperature, the interested reader may consult Ref. jljj]. For a general review of field 
theory at finite temperature, see Refs. 78- 8p(|. 



3.1 Evaluation up to order p 10 

In finite-temperature field theory the partition function is represented as a Euclidean 
functional integral 

Tr [exp(-H/T)] = J [dU] exp ( - J d A xC eff ^j . (3.1) 

The integration is performed over all field configurations which are periodic in the 
Euclidean time direction, U(x, x 4 + 0) — U{x,x^ with (3 = 1/T. The periodicity 
condition imposed on the magnon fields also reflects itself in the thermal propagator 

oo 

G(x) = A(z, 2* + 7l£) , (3.2) 



n = — oo 



where A(x) is the Euclidean propagator at zero temperature, 

frf J = 6(x 4 ) / /4 e^ k2x ^ Hx ^ . (3.3) 

(2tt) 4 ^-ih + fxH K } J (2tt) 3 1 J 
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An explicit representation for the thermal propagator, dimensionally regularized in 
the spatial dimension d s , is 

= Wf Q 2 ^ "T exp ( ~ 47^ ~ ^ Xn ) e(Xn) ' (3 ' 4) 

with 

x n = £4 + n/3 . (3.5) 

We restrict ourselves to the infinite volume limit and evaluate the free energy density 
z, defined by 

z = —T lim L~ 3 ln[Trexp(--H/T)] . (3.6) 

L—too 



In the evaluation of the various Feynman diagrams, we will repeatedly be dealing 
with thermal propagators (and space derivatives thereof), which have to be evaluated 
at the origin. It is convenient to introduce the following notation, 



G{x) 



G 



j 2=0 



AG(x) 



J x=0 



A n G{x) 



x=0 



(3.7) 



where A represents the Laplace operator in the spatial dimensions - no confusion 
should occur with A(x), which denotes the zero-temperature propagator. 

The quantities Gi, G&, as well as thermal propagators involving higher-order 
space derivatives, are split into a finite piece, which is temperature dependent, and a 
divergent piece, which is temperature independent, 



Gl = G T X + G[ 



G A 



G a + G , 



(3.8) 



The explicit expressions can be found in Ref. [36J and will not be given here. Rather, 
we would like to point out two important observations, which lead to a substantial 
simplification of the renormalization procedure. First, the temperature-independent 
pieces G^, G^, . . . are all related to momentum integrals of the form 



d ds k ( k 2 ) exp 



m 



0,1,2, 



which are proportional to 



exp[-x 4 iiH) / d s 
-. — 1 [m H 

( 7 x 4 ) m+ ^ V 2 



(3.9) 



(3.10) 



In dimensional regularization these expressions vanish altogether: G®,G\, and zero- 
temperature propagators involving higher-order space derivatives do not contribute 
in the limit d s — >3. 

The second observation concerns the fact that, up to order p 10 , the individual con- 
tributions to the free energy density from the various diagrams factorize into products 
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of thermal propagators (involving space derivatives or derivatives with respect to the 
magnetic field), which all have to be evaluated at the origin. As an example consider 
the two-loop graph 10a which yields the contribution 

zioa = -g|a (8*1 + 6/2 + 5/3) G A G A - ^ G x G A *. (3.11) 

According to the first observation regarding dimensional regularization, it is then clear 
that in the two products of thermal propagators above only the fully temperature- 
dependent pieces - G^^a an d Gj G\\ - are nonzero, whereas any other terms in- 
volving temperature-independent pieces of propagators vanish identically. We thus 
conclude that, using dimensional regularization, the renormalization of the partition 
function up to order p 10 is quite trivial. As we will see in the next subsection, the 
renormalization at the three-loop level, on the other hand, is more complicated, but 
still perfectly feasible within the effective field theory framework. 



Without going into more details (the interested reader may consult Ref. 36]), we 

present the final result for the free energy density of the 0(3) ferromagnet up to order 
p io. 



5 ~ e -»HnP 



3 2 2 
5^272 



105 / 9Z| 



3(8/! + 6/ 2 + 5/ 3 
128vr 3 S 2 7 5 




(3.12) 



The first term in this series does not depend on temperature and originates from the 
tree graph 2 (see Fig. [1]). The terms which involve half- integer powers of the temper- 
ature - T 5 / 2 , T 7 / 2 and T 9//2 , respectively - arise from the one-loop graphs displayed 
in Fig. [TJ They all contribute to the free energy density of noninteracting magnons. 
Remarkably, up to order p 10 , there is only one term in the above series - the contri- 
bution of order T 5 coming from the two-loop graphs 10a and 10b - which is due to 
the magnon-magnon interaction. 

In particular, there is no term of order T 4 in the above series for the free energy 
density: The two-loop graph 8, which would be the only candidate to yield such a 
contribution, is proportional to single space derivatives of the thermal propagator 
evaluated at the origin: 

z 8 <x[d r G(x)] x=0 [d r G(x)] x=0 = 0. (3.13) 

This contribution vanishes due to space rotation invariance of the leading-order effec- 
tive Lagrangian. 
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3.2 Evaluation at order p 11 



According to Fig. [5] we have a total of five diagrams at order p 11 . We first consider 
the two one-loop graphs which involve vertices from Cz f f, L\ n and C? e tt- For graph 



lid we obtain 



z \\d 



2di 



Ga 4 , 



yielding the temperature-dependent contribution 

rp 945 di 



647T2E7 2 



(3.14) 



(3.15) 



n=l 



n 2 



Graph lie is proportional to an integral over the torus T = TZ ds x 5 1 , with circle S 1 
defined by —(3/2 < x^ < (3/2, and involves a product of two thermal propagators, 



Zlle 



4/3C1 



d ds+l xA 2 G(x) A :i G(-x) 



(3.16) 



r 



This integral, however, can be reduced to an expression involving one propagator only, 



using the relation [36 



^(m+n) 



We then end up with 



ggpc) 
d(fiH) 



x=0 



J d ds+1 y A m G{-y) A n G{y) . (3.17) 



4/3C1 



A 



5 9G(x) 
d(fiH) 



x=0 



Accordingly, the temperature-dependent part of graph lie reads 



'lie 



10395 l 3Cl T n 



n=l 



n 2 



(3.18) 



(3.19) 



We now turn to the three-loop graphs - note that they exclusively contain vertices 
from the leading-order Lagrangian C 2 e ff- Graph 11a factorizes into a product of three 
thermal propagators (and space derivatives thereof), to be evaluated at the origin, 

zu a = -^^(d) 2 . (3.20) 

The subsequent three-loop graph lib, remarkably, does not contribute to the partition 
function, 

zii6 = 0. (3.21) 

As it was the case for the two- loop graph 8, the three- loop graph lib is identically 
zero. 
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Finally, for the cateye graph 11c we get 



zuc = +^G A (G 1 f. (3.22) 

The expression J stands for the following integral over the torus involving a product 
of four thermal propagators 

J = y d ds+l x d r G d r G d s G d s G , (3.23) 

where we have used the notation 

G = G(x) , G = G(-x) . (3.24) 

Note that the second term in f)3.22p cancels the contribution from graph 11a, such 
that the overall contribution from the three-loop graphs is the one proportional to the 
integral J. Remarkably, unlike all other pieces in the free energy density up to order 
p 11 , this quantity is not just a product of thermal propagators (or derivatives thereof) 
to be evaluated at the origin. The remaining task will be the renormalization and the 
numerical evaluation of this integral which contains a total of four infinite sums. In 
the next section and in Appendix |A] we address this problem in detail. 

Leaving aside these technical issues for a moment, we note that the cateye graph 
of order p u will lead to a term of order T 11 / 2 in the free energy density, 

JocT^. (3.25) 

Hence the spin-wave interaction in the low-temperature series of the free energy den- 
sity - beyond Dyson's T 5 -term - already manifests itself at order T 11 / 2 . It is re- 
markable that this contribution is exclusively determined by the symmetries of the 
leading-order effective Lagrangian £e// which involves the two couplings S and F - 
the spin-wave interaction at this order is not affected by the anisotropies of the cubic 
lattice. 



4 Renormalization of the Cateye Graph 

Using dimensional regularization it was straightforward to extract the finite pieces in 
the partition function up to two-loop order p 10 . The renormalization of the three- loop 
graph 1 1c, on the other hand, is more involved. We will follow the procedure outlined 
in Ref . [81] , where the same graph was considered in the context of a Lorentz-invariant 
effective field theory. 

To analyze the integral 

j= Jd ds+1 xd r Gd r Gd s Gd s G 
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in the limit d s — >3, we split the thermal propagator into two pieces 

G{x) = G T {x) + A(x) . 



(4.1) 



The ultraviolet singularities are contained in the zero-temperature propagator A(x), 
whereas the temperature-dependent part G T (x) is finite as d s — > 3. Note that, if we 
restrict ourselves to the origin, we reproduce the first relation of Eq. (I3.8j) . 

Inserting the above decomposition into the integral J, we end up with nine terms 
that can be grouped into the following six classeso 



A : 


G T (x) G T (x) G T (-x) G T (-x) 






B : 


A(x) G T (x) G T (-x)G T (-x), 


G T (x) G T (x) 


A(-x) G T (—x 


C : 


A 2 {x)G T (-x)G T {-x) , G T 


» G T (x) A 2 (- 


-x), 


D : 


A{x) G T {x)A{-x) G T {-x), 






E : 


A 2 (x)A(-x)G T (-x), A(x) 


G T (x) A 2 (-x) 




F : 


A 2 (x) A 2 (-x). 







(4.2) 

Terms of the classes D,E and F vanish identically since the product A(x)A(—x) of 
zero-temperature propagators involves the combination G(x4)G(— x^). The maximum 
number of B-functions a given term can contain - in order not to be zero - is two. 
Moreover, the arguments of the two G-functions have to coincide as it is the case with 
the terms of class C. We thus have to consider the cases A, B and C. 

The integral over the torus involving contributions of classes A and B, 
J d d ° +1 x (d r G T d r G T d s G T d s G T + 4 d r A d r G T d a G T d s G T } , 
converges at d s — 3. 



(4.3) 



Terms of class C, however, do lead to an ultraviolet-divergent integral. Consider, 
e.g., the term 

d r A{x) d r A{x) d s G T {-x) d s G T {-x) , (4.4) 

where we now have displayed the derivatives. For the zero-temperature piece d T A(x) 
we have 



x' 



d r A(x) (x — g- exp 

£42 



x 



47x4- 



(4.5) 



The Taylor series of the function d s G T (—x), evaluated at the origin, starts with a 
term linear in x, 

d s G T (-x) = d as G T (-x)\ x=0 x a + 0(x 3 ) . (4.6) 



2 For simplicity we do not display the derivatives. 
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Inserting this term into Eq.f l4.4p . we end up with the following contribution in J, 

J oc [d 3 xdxjAr) 2 e- £2/2 -< X4 x 2 , 



£4- 
f 1 

oc / dxi — 3- , (4.7) 

which is singular in the ultraviolet. On the other hand, one readily checks that this 
term in fact is the only one that has to be subtracted: The cubic Taylor term in the 
expansion of d s G T , Eq. fl4.6p . already leads to a convergent contribution to the integral 
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J. We now discuss the renormalization procedure in detail, along the lines of Ref. 
which we adapt to nonrelativistic kinematics. 

We first cut out a sphere S of radius |«S| < (3/2 around the origin and decompose 
the integral involving the contributions of class C according to 

d ds+l xd r Ad r Ad s G T d s G T 

= [d ds+1 xd r Ad r Ad s G T d s G T + [ d ds+1 x d r Ad r A d s G T d s G T . (4.8) 
Js Jr\s 

The integral over the complement T \ S of the sphere is not singular in the limit 
d s —t-3. In the integral over the sphere, which is divergent, we subtract the singular 
term discussed above, arriving at 

Jd da+l xd r A{x)d r A{x)d s G T {-x)d s G T {-x) 

= Jd ds+1 xd r A{x)d r A{x)Q ss {x) 

+Jd ds+1 * d r A(x)d r A(x) d as G T (-x) \ x=0 d Ps G T {-x)\ x=Q x a x? , (4.9) 

where the quantity Q ss (x) is defined as 

Q ss (x) =d s G T {-x)d s G T {-x)-d as G T {-x)\ x=Q d Ps G T {-x)\ x=0 x a xP . (4.10) 

Whereas in Eq. fl4.9p the first integral on the right hand side now is convergent, the 
second integral does contain the ultraviolet singularity. The last step in the isolation 
of this singularity consists in decomposing the respective integral as follows: 

J d d ° +1 x d r A(x) d r A(x) d as G T (-x) \ x=0 d Ps G T {-x) \ x=0 x a x? 

= [ d d ° +1 x d r A(x) d r A(x) d as G T (-x)\ x=0 d Ps G T (-x)\ x=Q x a x? 
Jn 

-! d d ^ 1 xd r A(x)d r A(x)d as G T (-x)\ x=0 d f}s G T (-x)\ x=0 x a x 13 . (4.11) 
Jn\s 
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The UV-singularity is contained in the integral over all Euclidean space, which can 
be cast into the form 

/ d d ° +1 xd r A(x) d r A{x) d as G T (-x)\ x=0 dp s G T (-x)\ x=0 x a x? 
Jn 



2 



4(4 + 2) _ rf , 2/ _ . d3 -2 e -»HnP d 

= zzi ' m L±* W 2 E r ( x - y) • ( 4 - 12 ) 

2 d<t « +t, 7r 2^/2 n=l n 2 j ^ 

In the limit 4 — > 3 the above regularized expression is finite and takes the value 

15 u^^l± e -^}\ (4.13) 



8192vr 4 7- \ ^ n% 

where we have defined the dimensionless quantity o as 



a = ^Hf3 = ^. (4.14) 



Collecting the various contributions we arrive at the following representation for 
the renormalized integral J: 

J = Jd A x (d r G T d r G T d s G T d s G T + 4 d r A d r G T d s G T d s G T 

+ 2 f d 4 x d r A d r A d s G T d s G T + 2 [d 4 x d r A d r A Q ss 
Jt\s Js 



2 f 4x9 r A5 r A9 QS G T (-x)U =0 «9^G T (-x)U = 
Jn\s 



a 

:Q X X^ 



'n\s 

15 „ii r- ^ e"™ 



l3 a/o- < — — i . (4.15) 
40967r 4 7 f V \~ n* J 

Note that all terms therein are well-defined at the physical dimension 4 = 3. 

Since the various integrands only depend on the variables r = \x\ and t — x^, the 
integrals become in fact two-dimensional, 

d 4 x = Anr 2 dr dt , (4.16) 

and the numerical evaluation of the integral J is straightforward. A very welcome 
consistency check on the numerics is provided by the fact that the result must be in- 
dependent of the radius of the sphere S. While more details concerning the numerical 
evaluation can be found in Appendix |AJ in the next section we discuss the result for 
the function J = J(cr). In particular, we consider the limit a— >0, which is needed for 
the evaluation of the spontaneous magnetization. 
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a 

Figure 3: The function j'(cr), where a is the dimensionless parameter a = fiH/T. 

5 Thermodynamics of the Ideal Ferromagnet 

For dimensional reasons, the renormalized integral J can be written as 

J{a)=T2 — ^ , <r = — i 7 = — , (5.1) 

>y 2 -t 2j 

where the quantity j'(cr) is a dimensionless function. A graph is provided in Fig. |3J 
In the limit a— >0, the function can be parametrized by 

3{°)=3i+32O + 0{a i ' 2 ). (5.2) 

The coefficients ji and j2 are pure numbers given by 

jx = 1.07 x 10" 5 , j 2 = -8 x 10" 5 . (5.3) 

It should be noted that, in the limit a — > 0, the last two contributions in Eq. (14.151) 
contain terms involving the square root yfa. Since they have opposite signs, however, 
they cancel each other. 

With the above representation for the quantity j(cr), the final result for the low- 
temperature expansion of free energy density for the ideal ferromagnet up to order 
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p u takes the form 



z 



0^272 n=1 n2 1D7T2L72 n=1 n2 

105 ( 9/ 2 \ . ^ e -^n/3 



^ ' \ ' / n=l 

3(8/! + 6/ 2 + 5/ 3 ) 5 ( A e-^' 
128tt 3 S 2 7 5 \h 

945 dj _n^e"^ n/3 IO395/3C1 11 ^ e"^ n/: " 
- — T~ \ 1 _ T~ > 

~A 3^ II / y I! ~ „, 3^o 13 ± / j 11 

64^2^7 2 ~ JJ 2 647T2L^ 7 2 ~ 72 2 



1 



2S 2 7 I 



+ C(T 6 ) . (5.4) 



Because the system is homogeneous, the pressure can be obtained from the temperature- 
dependent part of the free energy density, 

P = zq-z. (5.5) 

Accordingly, up to order p 11 , the low-temperature series for the pressure reads 

P = h T 5 * + h x T* + h 2 T^ + h 3 T 5 + h 4 T% +C(T 6 ), (5.6) 

where the coefficients hi are given by 

k = " 3 3 2^ 1 ' 

8^272 n=1 n2 

15 / 3 ^ e -^n/3 
"' 1 ~~ 77, 3 7 2^ 2 ' 

107T2L72 «2 

105 / 9/ 3 2 \ ~ 
2 " 327rls 7 I ^£7 C 7 ^ „i ' 

3(8/i + 6/ 2 + 5/3) f A e-» Hn P 



h« = 



128tt3S 2 7 5 {5 nf 



945 / lU 3 d\ ^ e-^ 1 . 

^4 = 3 ^ 11 \di- -=— 2^ — ~ + ^r^J- ( 5 - 7 ) 

647T2S7^" y ^7 J ~l TL 2 ^ 

In the limit a = jiH/T— >0, these coefficients become temperature independent and 
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the sums reduce to Riemann zeta functions, 



1 



"'! - ± 3 C(s) ' 

8^272 

I - ^3 A ,J, 

hl ~ ~ 3 7 (,{2) 

1D7T2 L72 



105 / 9/f . 



3(8li + 6Z 2 + 5/3) 2 5 , 
128tt 3 E 2 7 5 C l2 ' 

945 f 4 _nta, {( „ )+ - (5 . 



647r2E 7 - \ £7 / 2E 2 7 t 

The spin-wave interaction manifests itself in the last two terms involving the coeffi- 
cients /13 and /14. The contribution proportional to five powers of the temperature in 
the pressure is the famous Dyson term. In the effective theory it originates from the 
two-loop graphs 10a and 10b of Fig. [TJ Our main new result concerns the manifesta- 
tion of the spin- wave interaction beyond Dyson: the leading correction in the pressure 
is of order T 11 / 2 . It is contained in the last term of the coefficient J14 and stems from 
the three-loop graph 11c. 

Note that all other contributions in the pressure up to order p 11 originate from 
one-loop graphs - those graphs describe noninteracting magnons and merely modify 
the dispersion relation. In the above series for the pressure, they involve half-integer 
powers of the temperature: T 5 / 2 , T 7 / 2 , T 9 / 2 and T 11 / 2 . 

We have to point out that the sign of the Dyson term of order T 5 is not determined 
by the symmetries - the low-energy constants l\, I2 and I3 appearing in the coefficient 
h 3 may take positive or negative values, depending on the specific underlying model. 
For the present case of the Heisenberg model, however, Dyson has derived an explicit 
microscopic expression for /13. As it turns out, for all three types of cubic lattices, 
this coefficient is positive, leading to a positive contribution to the pressure. We thus 
conclude that the spin-wave interaction in the ideal ferromagnet is repulsive at low 
temperatures. 

Remarkably, while the sign of the coefficient of order T 5 is not determined within 
the effective theory framework, the sign of the coefficient of the subsequent interaction 
contribution of order T 11 / 2 is unambiguously fixed: the last term in only involves the 
coupling constants of the leading-order effective Lagrangian C 2 e ^ and the coefficient j\ 
which is a pure number - the conditions imposed by symmetry are thus very restrictive 
here. Since the numerical value of ji is positive, the corresponding contribution to 
the pressure is positive as well, enhancing thus the weak repulsive interaction between 
spin waves at low temperatures. 
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Finally, let us consider the low-temperature series for the energy density u, for the 
entropy density s, and for the heat capacity cy of the 0(3) ferromagnet. They are 
readily worked out from the thermodynamic relations 

In the limit a— >0, we obtain 



u = ~h T% + \~h x Tl + fh 2 T% + 4h 3 T 5 + fh 4 T^ + £>(T 6 ), 



s = lh T2 + \h x T* + \h 2 T* + 5/i 3 T 4 + f /j 4 T2 + C(T 5 ) , 
c y = f fr T§ + f ^Tl + f h 2 T 7 * + 20/* 3 T 4 + f h 4 Tl + 0(T 5 ). (5.10) 



Again, the correction to Dyson's result is contained in the respective last terms in the 
above series involving the coefficient h 4 . 



6 Spontaneous Magnetization: Effective Framework 
versus Condensed Matter Literature 



We now turn to the discussion of the general structure of the low-temperature series 
for the spontaneous magnetization of the ideal ferromagnet. While this problem has 
attracted more than a hundred authors over the last few decades, to the best of our 
knowledge, a rigorous and fully systematic calculation of higher-order corrections to 
the Dyson term has never been achieved. Before we review the relevant results in 
the literature, let us analyze the problem within the systematic effective field theory 
framework. 

With the expression for the free energy density (15. 4p . the low-temperature expan- 
sion for the spontaneous magnetization 

(~)7 

E(T) = - lim — — (6.1) 
h^o d(fiH) 

of the 0(3) ferromagnet, up to order T 9//2 , takes the form 

= 1 - a T% - a x Ti - a 2 T^ - a 3 T 4 - a 4 T% + C(T 5 ) . (6.2) 
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The coefficients cti are independent of the temperature and given by 



"0 = ~ 3^ 3 C(|), 

57T2 L72 
107T2 2V72 

105 / 9Z§ \ A/7x 

3(8/! + 6/ 2 + 5/3) 5 3 

«3 " 647r 3 E 3 7 5 CCaKCa). 

945 / ll/ 3Cl \ 9 1 . 

"4 = 3 U Ni - — — C(| - r7^J2 ■ (6.3) 

647r2E 2 7~ y ^7 / 2S 3 72 

Up to order T 4 , we reproduce Dyson's series. In the effective Lagrangian framework, 
the famous interaction term of order T 4 in the spontaneous magnetization originates 
from the two-loop graphs 10a and 10b which involve vertices from the next-to-leading 
order Lagrangian £ 4 ^. Note that there is no interaction term of order T 3 in the 
above series - for many years such a spurious term has haunted the condensed matter 
literatureH 

Our main new result concerns the leading correction to Dyson's term, which origi- 
nates from the three-loop graph 11c. The correction in the spontaneous magnetization 
is of order T 9 / 2 . Remarkably, the corresponding coefficient - the last term in - 
does not involve any higher-order low-energy constants. It only involves S and F, as 
well as the coefficient j'2, which is a pure number determined by the symmetries of the 
underlying Heisenberg model. Since the coefficient 32 is negative, this contribution has 
the same sign as the Dyson coefficient 03. The effect of the three-loop contribution is 
thus to enhance the weak spin-wave interaction found by Dyson. 

Apart from these two interaction terms of order T 4 and T 9//2 , respectively, all 
other temperature-dependent contributions to the spontaneous magnetization orig- 
inate from one-loop graphs, which describe noninteracting magnons. They merely 
modify the dispersion relation or - as Dyson expressed it [6( - they merely arise from 
the discreteness of the lattice, are easy to calculate and are not of any theoretical in- 
terest. In the above series for the spontaneous magnetization, they involve half-integer 
powers of the temperature: T 3/2 , T 5/2 , T 7/2 and T 9/2 . 

It should be pointed out that the contribution of order T 9 / 2 contains two parts: 
The first term in the coefficient «4 is due to two one-loop graphs associated with non- 
interacting magnons. The second term is due to a three-loop graph and represents the 



3 The first encounter with this spurious term seems to date back to the year 1958 (see Ref. [19(). 
It then continued haunting the literature over a period of at least 25 years until 1983 when it was 
last sighted in Ref. [29| . 
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dominant spin-wave interaction term beyond Dyson. Note that the Dyson coefficient 
0:3, on the other hand, exclusively involves an interaction part. 



Here comes the appropriate place to compare our results with the condensed mat- 
ter literature. Indeed, several authors - most notably, Dyson himself - also have 
discussed the structure of the low-temperature series for the spontaneous magnetiza- 
tion beyond order T 4 . We make our comparison along four lines of observations. 

Our first observation is that all published calculations or estimates of higher-order 



interaction terms [q. I82H84] apparently failed to identify the dominant T 9//2 -correction 



to the Dyson term in the spontaneous magnetization. 

The second observation is that there appears to be consensus in the literature on 
how graphs related to the iwo-spin-wave problem should manifest themselves beyond 
T 4 . Dyson classified his terms according to the quantity F, where F is the number 
of independent particles which are concerned in the interactions which the particular 
term describes |6|]. For F = 2, which is referred to as the two-spin-wave problem in 



Ref. [821 ]. the corresponding corrections are expected to show up at order T 5 according 



to Refs. 82- 84] - hence these authors seem to agree on that the dominant correction 
to the Dyson term should be of order T 5 in the spontaneous magnetization. However, 
this claim is not correct - it is in contradiction with the fully systematic effective field 
theory analysis which has demonstrated that the dominant correction sets in at order 
T 9 / 2 . 

The third observation concerns the i/jree-spin-wave problem, i.e. the effect of in- 
teraction terms with F — 3. Dyson identified two such contributions - formulae (128) 
and (130) in his second article of Ref. 0] - and showed that in the spontaneous mag- 
netization these are of order T 13//2 and T 5 , respectively. In the article by Morita and 
Tanaka [82| , however, it is claimed that the i/jree-spin-wave problem starts manifesting 
itself at order T 13 / 2 , missine; thus the term of order T . 

Finally, the fourth observation is that the only place in the literature where a Feyn- 
man diagram displaying the cateye structure of graph 11c seems to have appeared, is 



in the more recent article by Chang 83]. However, he concludes that interactions orig- 
inating from such a diagram start showing up only at order T 15//2 in the spontaneous 
magnetization. This claim, again, is erroneous, as it contradicts the fully systematic 
effective theory analysis which has demonstrated that the leading term originating 
from a cateye graph is of order T 9 / 2 . 

In view of the quite impressive collection of temperature powers established over 
the years, one may easily get confused - after all, may some of these temperature 
powers, again, merely be spurious? One would certainly like to gain some deeper 
insight into the general structure of the low-temperature series for the spontaneous 
magnetization beyond the leading correction to the Dyson term. Let us therefore ad- 
dress the problem in a fully systematic way within the effective field theory framework 
- first on the level of the free energy density. 
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12a 



13a 



Figure 4: Two Feynman graphs related to the low-temperature expansion of the 
partition function for a ferromagnet at order p 12 and p 13 in dimension d—3+1. The 
numbers attached to the vertices refer to the piece of the effective Lagrangian they 
come from. Note that there are further Feynman graphs of order p 12 and p 13 which 
we have not displayed. 

Indeed, it is quite easy to see that corrections due to the spin-wave interaction 
continue to proceed in steps of T 1//2 . In Fig. H] we have displayed some of the relevant 
higher-order graphs which contribute beyond order p 11 or, equivalently, beyond T 11 / 2 . 
At order T 6 in the free energy density, the two-loop graph 12a with an insertion from 
£e// contributes, while at order T 13 ^ 2 the three-loop graph 13a with a vertex from 
C^jf is relevant. 

We may classify the graphs of the effective theory according to the number of 
loops they contain and discuss the various contributions at the level of the spon- 
taneous magnetization. Interactions related to two-loop diagrams start manifesting 
themselves through the Dyson term of order T 4 and then proceed in integer steps of 
T - these are the two- loop graphs with successive insertions from ^e//> ^e//> • • • > 
giving rise to terms of order T 4 , T 5 , T 6 , ... in the spontaneous magnetization. Inter- 
actions related to three-loop diagrams start showing up at order T 9//2 . They give rise 
to the dominant correction to Dyson's result, and then also proceed in steps of T - 
in the effective theory these correspond to the three-loop graphs with successive in- 
sertions from Clff, £>tff, £tffi ■ ■ ■ i leading to half-integer powers of the temperature: 
T 9 / 2 , T 11 / 2 , T 13 / 2 , . . . Four- loop interactions are expected to enter the game at order 
T 6 . They will continue contributing to the spontaneous magnetization in ascending 
powers of T through terms of order T 6 , T 7 , T 8 , . . . 

Accordingly, the low-temperature expansion for the spontaneous magnetization of 
the ideal ferromagnet exhibits the following general structure: 

= 1 - a T% - ai ri - a 2 T 7 2 - a 3 T 4 - ct 4 Ti - ct 5 T 5 - a 6 T^ + C(T 6 ) . (6.4) 

Note that we have highlighted all terms which are related to the spin-wave interaction. 
We thus realize that the various published temperature powers as such are not in 
contradiction with the effective field theory prediction - sooner or later, the effective 
expansion will hit them all. The point is, however, that there are two gaps in the 
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multitude of published temperature powers - interaction terms of order T 9 ' 2 and 
T 11 / 2 , to the best of our knowledge, have never been identified. 

One may say that the effect of the spin-wave interaction on the low-temperature 
series for the spontaneous magnetization of an ideal ferromagnet is quite peculiar. On 
the one hand, the interaction starts manifesting itself only at order T 4 , i.e. far beyond 
the Bloch term of order T 3 / 2 . On the other hand, subsequent interaction corrections 
closely follow the term of order T 4 , starting with T 9//2 and then proceeding in steps 
ofT 1 / 2 . 



7 Conclusions 



The question of how the spin-wave interaction manifests itself in the low-temperature 
expansion of the spontaneous magnetization of an ideal ferromagnet, has a long his- 
tory. Early attempts that ended up with temperature powers of order T 7 / 4 and T 2 
turned out to be wrong altogether, as shown by Dyson's rigorous analysis which 
demonstrated that the spin- wave interaction sets in only at order T 4 . 

After Dyson's rather complicated analysis, there emerged an active phase of re- 
search in which many authors tried to derive the T 4 -term in the spontaneous mag- 
netization in a more accessible manner. It is interesting to note that this problem 

- using conventional condensed matter methods - indeed seems to be nontrivial, as 
various authors all of a sudden ended up with an interaction term of order T 3 in 
the spontaneous magnetization, thereby contradicting Dyson. As it turned out, these 
attempts were wrong altogether, as they were plagued with errors originating from 
the approximate methods used. 

We would like to stress that the fully systematic effective Lagrangian method, 
based on symmetry considerations, does not display any such defects. Rather, the 
existence of Dyson's T 4 -term - and, at the same time, the absence of a T 3 -term 

- is an immediate consequence of the underlying symmetries inherent in the ideal 
ferromagnet. 



In his collection of selected papers that appeared in 1996, Dyson comments [181 ] : 
"After 1966, the subject of spin-wave interactions went into a long sleep. "0 Still, so 
it seems, sporadically over time, the peaceful sleep has been interrupted, as several 
authors were attracted by the problem of how the general structure of the series 
for the spontaneous magnetization of an ideal ferromagnet beyond the Dyson term 
should look like. Remarkably, not all of the various findings are consistent with one 
another. It was our motivation to solve this paradox, making use of the systematic 
and model-independent method of effective Lagrangians. 



*In 1966 Keffer's comprehensive review on spin waves appeared (Ref. [33[). 
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As we have demonstrated in the present study, it is rather straightforward to 
go beyond Dyson's analysis by considering three-loop effects in the effective field 
theory. Still, the explicit evaluation at order p u is quite non-trivial as it involves the 
renormalization and subsequent numerical evaluation of a three-loop graph, which is 
proportional to an integral over a product of four thermal propagators, each one of 
them involving an infinite sum. The corresponding interaction term beyond Dyson is 
completely fixed by the symmetries of the leading-order effective Lagrangian C 2 e ^ - 
lattice anisotropies, showing up at higher orders in the effective Lagrangian, do not 
affect this result. What is quite remarkable is the fact that all previous attempts to 
go beyond Dyson apparently have failed to correctly identify this interaction term of 
order T 9 / 2 in the spontaneous magnetization of an ideal ferromagnet. 

We have also discussed the origin and the structure of even higher-order corrections 
in the low-temperature expansion of the spontaneous magnetization, pointing out that 
they continue to proceed in steps of T 1//2 beyond the contribution of order T 9//2 . Again, 
earlier attempts to gain insight into the general structure of this series were incorrect. 

The present study has thus solved - once and for all - the problem of how the spin- 
wave interaction in an ideal ferromagnet manifests itself in low-temperature expansion 
of the spontaneous magnetization beyond the Dyson term. 

Hopefully, we have convinced the reader that the effective Lagrangian technique 
does not merely consist in rederiving known results or in rephrasing condensed matter 
problems in another language - rather, in many cases as in the present one, it clearly 
proves to be more powerful than conventional condensed matter methods, allowing 
one to go to higher orders of the low-temperature expansion in a controlled and sys- 
tematic manner. In view of the many articles that have dealt with the problem of the 
manifestation of the spin-wave interaction in an ideal ferromagnet at low temperatures 
- concerning both the extension of Dyson's series to higher orders and the confusion 
regarding the spurious T 3 -term - it is quite striking how efficiently the effective theory 
analysis settles all these questions in a conclusive way. 

We do not claim to have contributed to the actual experimental situation - spin- 
wave interactions in a ferromagnet are very weak. Above all, there are many interac- 
tions in addition to the exchange interaction in a real ferromagnet, which would also 
have to be accounted for in a more realistic approach. While this would perfectly be 
feasible within the effective Lagrangian framework, here we have restricted ourselves 
to the simple model of the ideal ferromagnet - after all, it is for this system where 
corrections to Dyson's result have been derived over the years. 

Although there is no such object as a perfectly ideal ferromagnet in nature, still, 
the "clean" ideal ferromagnet could be investigated in a numerical simulation of the 
Heisenberg model and the existence of the T 9//2 -term in the low-temperature expansion 
of the spontaneous magnetization might be verified this way. 
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A Numerical Evaluation of the Cateye Graph 

To numerically evaluate the integral J defined in Eq. ( 14. 15ft . we introduce the dimen- 
sionless variables rj and £, 

i [r 

r] = Tx 4 , £=2V ^ A '^ 

In the integrals over the torus which involve quartic and triple sums - the first two 
terms in Eq. (I4.15P - we first integrate over all three-dimensional space, ending up 
with one- dimensional integrals in the variable rj. For the quartic sum we obtain 



d 4 x d r G T (x) d r G T (x) d s G T (-x) d s G T {-x 

is ii r 1 ' 2 



r 

T 



2048vr9/ 2 7 1 3/2 J ar ' 2^ 

1 ^ -1/2 ni...n 4 =l 



-7/2 



1 



ri+m rj+n 2 —r)+n 3 — ri+714 

Q(v,ni,n2,n 3 ,n 4 ) = — — , (A.2) 

Mr) + ni)(i) + n 2 )(-r) + n 3 ){-r) + n 4 )J 

while for the triple sum we get 

d 4 x d r A(x) d r G T {x) d s G T {-x) d s G T {-x) 

15 « r 1/2 



20487T 9 / 2 7 13 / 2 



TT / dr] e ~ a 

^° n 2 -riA=l 



(n 2+ n 3+ n 4 ) Q^, , n 2) n 3 , n 4 ) 

ri2...ri4 = 

7/2 



7; 17+712 — ?7+7T.3 -r7+n 4 

Q(?7, 0, n 2 , n 3 , n 4 ) = 7 ; 5/2 , (A.3) 



^ + n 2)(~V + n 3 )(-r) + n 4 )J 



with 
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Note that for the triple sums the integration over rj only extends over the interval 
[0, |], due to the G-function contained in the zero-temperature propagator A(x). 

The quantities Q(r], ni, n 2 , n 3 , 714) and Q(rj, 0, n 2 , n 3 , n 4 ) depend in a nontrivial 
manner on the summation variables. The slowest convergence for the entire expres- 
sions Eq. (IA.2j) and Eq. flA.3j) is observed for the case a = 0, where no exponential 



damping occurs. We have performed the numerical summation in a " Cartesian" way. 
We first define the vector iVj = (n 1; n 2 , n 2 , n 4 ). The first partial sum Si in the quartic 
series simply corresponds to the combination N\ = (1, 1, 1, 1) of indices. The second 
partial sum S 2 then contains all combinations of indices in the vector N 2 with at least 
one index equal to two: (2, 1, 1, 1), ... , (2, 2, 2, 2), etc. For large values of % and for 
a = 0, the partial sums Si converge according to 1/ S^ 2 . Proceeding in an analogous 
manner for the triple sums, one obtains the same asymptotic behavior. 

Expressions suitable for the numerical evaluation of the remaining three integrals 
of Eq. (I4.15P involving double sums are 

d 4 x d r A{x) d r A{x) d s G T {-x) d s G T {-x) 

T\S 

S poo 



1 po poo _ _ 

19a 5 13/2 t¥ / dv I dtt e J2 e^ (ni+n2) P(t 



ri,n x ,n 2 ) 



_^2| 2_ L [ 



7] — 77+711 — 71 + 712 



P(£,77,ni,n 2 ) = (A.5) 

T] 2 (-r] + ni)(—r} + n 2 ) 



<d 4 xd r A(x)d r A(x) Q ss (x) 



(A.6) 



1287r 5 7 i3/2 T 2 J o dr l ju 



V^ 1 



d^e J2 e-^ +n ^ Q(lr, V , ni ,n 2 ,a), 

ni ,ri2=l 



with 



Q(C,V,ni,n 2 ,a) 





+ 1 4 


- i 


e 2 ^ e ?2( 


i , i n 


— 7] + ni 


-17 + 772 I 


— Tj + 7l-\ —Tj + 712 






[{(-V+ni)(-V+n 2 )} 5/2 


5/2 5/2 
V n 2 -1 



(A.7) 
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and finally, 



d 4 x d r A(x) d r A(x) d sa G T (-x)\ x=0 x a d sP G T (-x) \ x=0 

■y poo poo 



poo poo ^ 

Js Jo 



128^7-/2 J s JQ 



e 



-2e/v 



R(£, 77, m, n 2 ) = -5/2 • ( A -8) 

It is understood that in the above integrals the radius of the sphere is chosen as S = \. 
For large values of i and for o = 0, the partial sums Si related to the above three 
expressions involving double sums also converge according to 1/S^ 2 . 
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